Localized Perturbations of Integrable Systems 
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The statistics of energy levels of a rectangular billiard, that is perturbed by a strong localized 
potential, are studied analytically and numerically, when this perturbation is at the center or at a 
typical position. Different results are found for these two types of positions. If the scatterer is at 
the center, the symmetry leads to additional contributions, some of them are related to the angular 
dependence of the potential. The limit of the 5-like scatterer is obtained explicitly. The form factor, 
that is the Fourier transform of the energy-energy correlation function, is calculated analytically, in 
the framework of the semiclassical geometrical theory of diffraction, and numerically. Contributions 
of classical orbits that are non diagonal are calculated and are found to be essential. 



PACS numbers: 03.65.Nk, 03.65.Sq, 05.45.Mt 

The distribution of energy levels exhibits a high de- 
gree of universality and is a central subject in the field 
of "Quantum Chaos" |^ . For systems that are chaotic 
in the classical limit the statistics are of Random Matrix 
Theory (RMT)) while for typical integrable systems 
the level distribution satisfies Poissonian statistics Q . In 
the semiclassical regime this universal behavior holds for 
a wide range in energy. There are also regimes of energy 
where spectral correlations related to periodic orbits are 
important 1^, ^. In intermediate situations such a high 
degree of universality in not found. For mixed systems, 
where in some parts of phase space the motion is chaotic 
and in other parts it is regular, the statistics exhibit some 
general features [0, || . Another type of intermediate be- 
havior may be found for integrable systems perturbed 
by singularities of spatial extension that is much smaller 
then the wavelength of the quantum particle. Examples 
of relevant systems are billiards with flux lines, sharp cor- 
ners and 5-like interactions |l^, |ll|. Here we report 
results obtained for a rectangular billiard perturbed by a 
(5-like impurity |l^ , known as the Seba billiard |^ . Some 
of these results can alternatively be concluded from a re- 
cent general formulation by Bogomolny and Giraud [ p^ . 

The interest in billiards of various types is primarily 
theoretical since it is relatively easy to analyze them an- 
alytically and numerically. Billiards were studied also 
experimentally for electrons , microwaves |]l5| and for 
laser cooled atoms Q. We hope that in the future, per- 
turbations of the type discussed in the present work will 
also be introduced experimentally. 

Trace formulas that express the quantum density of 
states in the semiclassical limit as sums over classical 
periodic orbits were derived for chaotic JT?] , and in- 
tegrable systems. For perturbations smaller than 
the wavelength, standard semiclassical theory used in 
the derivation of these formulas fails and diffraction ef- 
fects have to be taken into account. This can be done 
in the framework of the Geometrical Theory of Diffrac- 



tion (GTD) |20j. In this approximation, which is valid 
far from the perturbation, the Green's function for the 
system (without the boundary) is given by G{k; r, r') ~ 
Go(fc; r, v')+Go{k- r, Vo)D{e, e')Go{k; Tq, r'), where 6 and 
6' denote the directions of r — ro and rg — r' respectively 
and Go is the free Green's function. The diffraction con- 
stant, D, describes the scattering from the perturbation. 
For the rectangular billiard with a (5-likc perturbation, 
that is subject of the present work, only diffraction effects 
are responsible for the deviations from the behavior of in- 
tegrable systems. Therefore this is an ideal system for the 
exploration of such effects. Moreover for this problem the 
analytical and numerical calculations are relatively easy. 
The statistics depend on the location of the scatterer and 
on the boundary conditions 1 1 1 , pli . This is in contrast 



to chaotic systems where the spectral statistics are not 
affected by such scatterers . 

The diagonal approximation where only contribu- 
tions from orbits with equal actions are considered is ex- 
tensively used in the field of "Quantum Chaos". It is 
not applicable for systems with localized perturbations. 
A method to take into account dominant non diagonal 
contributions, in integrable systems, was developed by 
Bogomolny and will be used here. 

In this work a rectangular billiard with sides and ay, 
such that the aspect ratio is irrational, perturbed by a 
localized scatterer is studied. The scatterer is represented 
by a potential of typical size a such that 
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where /(y) is small where y is large. 

The diffraction constant is the on shell matrix element 
of the T matrix 0(9', 9) = (k|T(£')|q) where k is the out- 
going momentum (in direction 9') and q is the incoming 
momentum (in the direction 9). The energies of the in- 
coming and outgoing waves are equal, that \sk = q ~ \/E 
(in units ?i = 1 and to = ^ used in this letter). The Born 
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series cannot be used to compute D{9',9) when fca <C 1 
since the free Green's function diverges as In ka at short 
distances. A method that is regular when fca <C 1 was 
introduced by Noyce j2^. In this method the scattering 
in the forward direction is resummed p^ . It leads to a 
diffraction constant that is a ratio of series. The series 



in the numerator and in the denominator are expanded 
in the number of scattering events (just like the Born 
series). Every term in these series is then expanded for 
ka ^ 1 (up to terms of order fc^a^) and both series are 
summed (with respect the number of scattering events) 
to give the angle dependent diffraction constant 
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c.c. denotes the complex conjugate, 7 is Euler's constant 
and = / d2y/(y). Also V(l), 0(1), Mo, Ml, M2 and 
Avicd are constants, independent of 9, 9', and Mi, M2 
and Alcd depend logarithmically on ka. These constants, 
given by series of integrals, involving the potential (|^), 



were calculated in |12 

In the limit a and k fixed, a finite diffraction 
constant is obtained if the potential is such that, ^^^^ ~ 
^[\n{a/l) + B], where I and B are constants, leading to 

i^^2,r(^-ln(M^-7 + s)". (4) 

It depends on the combination _B — In / of the two param- 
eters / and B. Therefore these are somewhat arbitrary 
in the limit a 0. 

First we assume that a is sufficiently small so that (|^) 
can be approximated by (|^), and because of its slow varia- 
tion with fc, it can be replaced by the constant D. The os- 
cillatory part of the density of states, in the semiclassical 
limit, is a sum over contributions of periodic and diffract- 
ing orbits. Diffracting orbits are orbits which start and 
return to the scatterer. For the rectangular billiard with 
a localized (angle independent) scatterer at its center the 
density of states is |Q 
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■ , where xi ~ Nj^ + Alj^ , 

and Xi ~ Xi^i + Nj^ + Avij. . The area of the billiard is A, 
the length of a periodic orbit is Ip and Ij is the length of 
a diffracting segment j with Nj and Afj reflections from 
the boundary. The density of states (|^) that is expanded 
to the third order in D, is used to compute the correlation 
function 
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and its Fourier transform, the form factor 

/oo 
dr,R^{il)e^^^^\ (7) 
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where A denotes the mean level spacing. The brackets 
denote averaging over an energy scale much larger then 
A but much smaller then E. If only the contributions 
from periodic orbits to the density of states are taken 
into account the form factor is given by: 

i;j.^^,^,,,..-v.,(|._i±v)) 

(8) 

where Sp is the action of the orbit, tp is its period, 
and r > is assumed. The diagonal approximation 
can be used to compute the contributions from peri- 
odic and once diffracting orbits. When there are more 
then 3 segments of orbits in the exponent of (|^), non di- 
agonal contributions are of importance, since then one 
finds saddle manifolds consisting of different combina- 
tions of orbits with almost identical total length so that 
their phase is almost stationary An example of 

such a saddle manifold is given by a periodic orbit of 
length Ip — 2J Nptt'^ + MpUy, and the pairs of diffract- 



ing segments of lengths Ij^ = ^ N'j a^ + A'hj a^ that sat- 

isfy N^, + N^, = 2Np, M,, + M,, = 2Mp and ^ ~ ^. 
The length difference Ip — Ij^ — Ij^ is small, of the order of 
1/k. Since in billiards the action is Sj = klj the action 
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difference is of order unity and these contributions are in 
phase. Other non diagonal contributions of this type can 
contribute significantly as well. The resulting form factor 
for a scatterer at the center, up to order t^, is found to 
be 



K{t) = 1 
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To obtain @ we used the optical theorem, that for angle 
independent scattering is 
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If the scatterer is at a typical location, namely shifted 
from the center by {5xax,Syay), with 6x, 5y and S^/Sy all 
irrational, the form factor is [|l2| 
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The difference between (||) and ( |ll| ) is due to length de- 
generacies. When the scatterer is at the center there are 
four diffracting segments of identical length, while if it is 
moved from the center this degeneracy is broken. For a 
quarter of all diffracting segments, for which Nj and Mj 
are odd, the degeneracy is totally lifted. For orbits with 
even Nj and Mj the location of the scatterer does not 
affect this degeneracy. For the rest of the segments the 
degeneracy is only partly lifted. When all length degen- 
eracies are taken into account one obtains (pi]). 

The form factor can be compared with numerical re- 
sults obtained for the case of point interactions, where the 
eigenvalues are the roots of some function, and therefore 
can be easily found numerically . The form factor was 
calculated for several values of D and compared to the 
analytical result (|Tl|) in Fig. |l|. Agreement with (^l|) is 
found for short times, as is expected. 

For a scatterer at the center only levels with wave func- 
tions that are symmetric with respect to the X and Y 
axes are perturbed by the scatterer. Since the value 
of these wave functions for all eigenvalues is the same 
(i/)„(x = 0) = the resulting equation is the same 

as for the Seba billiard with periodic boundary condi- 
tions The form factor of the perturbed levels is 
related to the one of the full spectrum. The eigenvalues 
of the four different symmetry classes of the rectangle 
can be assumed to be uncorrelated, leading to 
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(4r). 



(12) 



The form factor calculated from all levels and the scaled 
form factor obtained from the perturbed levels with the 
help of (|l^) are compared to the analytical result (H), 
for D = -Ai, in Fig. ||. It is clear that (|l2|) is vahd. 
For very small times the full form factor deviates from 
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FIG. 1: The form factor of a scatterer at a typical position, 
for some diffraction constants, thin lines, compared to the 
analytical result heavy lines. 
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FIG. 2: The form factor for the scatterer at the center 
(xo = 0) compared to the scaled form factor calculated from 
perturbed levels (jl^) and the analytical result ^ 



its expected value since there are not enough orbits that 
contribute (the calculation is not semiclassical enough). 

In order to study the effect of the angle dependence of 
the diffraction on the form factor it was calculated to the 
order for the diffraction constant (|^). For a scatterer 
at a typical location the form factor is found to be [n2| 



K{t) = 1-\\C\^t 
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This form factor is similar to (pd]). Since C of (^) sat- 
isfies the optical theorem (p^), this form factor can also 
be obtained from an angle independent potential, with 
the diffraction constant C. If the scatterer is at the cen- 
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ter 1^ 



i^(r) = l-l|qV+i|C|W 



(14) 



with C" = 1 — 2fc^a^Mo, where Mq is related to integrals 
over the potential (0). It resembles the form factor (||) 
that was obtained for angle independent scattering. The 
modification is of the order k^a^ and typically cannot 
change the sign of the expansion coefficients. 

The condition for the applicability of the approxima- 
tions used in this letter is a <ti X ^ a^, CLy, where A 
is the wavelength of the particles . Up to corrections of 
order (ka)^, the form factor ( [l^ ) reduces to (|l|) in the 
order r^. Therefore the angle dependence plays no role 
up to this order. If the scatterer is at the center the sit- 
uation is somewhat different as can be seen comparing 
( |r^ ) with (H). There is a correction C resulting of the 
angular dependence of D{9' , 0) given by (||). It is a conse- 
quence of the increased number of length degeneracies of 
the diffracting orbits when the scatterer is at the center. 
Since the form factor (^3|) describes essentially angle in- 
dependent scattering the limit a ^ describes correctly 
the physics of the regime a < A. This is so although the 
classical dynamics (in the long time limit) are expected to 
be chaotic in nearly all of phase space and similar to the 
ones of the Sinai billiard. This robustness improves the 
chances for the experimental realization of the results of 
the present work. For a ^ A, semiclassical theory works 
and the system should behave as a Sinai billiard, with 
level statistics given by Random Matrix Theory (RMT) 
lH, H (with deviations, see j25|). 

The spectral statistics found in the present work dif- 
fer from the ones of the known universality classes. It is 
characterized by the form factor of the type presented in 
Figs. |l| and |[ This form factor is equal to 1 at t = 0, 
resulting of the fact that for small r the number of clas- 
sical orbits that are scattered is small. The contribution 
that is first order in r originates from the combinations 
of forward diffracting orbits and periodic orbits. These 
always have the same lengths leading to the contribution 
SDr. By the optical theorem ( p^ ) it is always negative. 
For T 3> 1 the form factor approaches unity because of 
the discreteness of the spectrum This general de- 
scription should hold for other integrable systems, that 
are perturbed similarly. 
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